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Abstract. It has been recently shown that the Bethe-Larkin formula for the energy 
losses of fast heavy-ion projectiles in dense hydrogen plasmas is corrected by the 
electron-ion correlations [Phys. Rev. Lett. 101, 075002 (2008)]. We report numerical 
estimates of this correction based on the values of <? e i(0) obtained by numerical 
simulations in [Phys. Rev. E 61, 3470 (2000)]. We also extend this result to the case of 
projectiles with dicluster charge distribution. We show that the experimental visibility 
of the electron-ion correlation correction is enhanced in the case of dicluster projectiles 
with randomly orientated charge centers. Although we consider here the hydrogen 
plasmas to make the effect physically more clear, the generalization to multispecies 
plasmas is straightforward. 
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1. Introduction 

The calculation of the stopping power for a projectile traversing a plasma can be 
performed within the so called polarizational approach. The Lindhard formula relates 
the polarizational stopping power with the (longitudinal) dielectric function of the 
medium pQ. When a projectile with a charge distribution p(r,t) moves at a constant 
velocity v through an isotropic Coulomb fluid, this formula can be written as [2] 

"f = ^/dk^(-Im £ -'ft k -v)) W k)P, (1, 

where k = |k|. The polarizational approach to the stopping power calculation is justified 
provided that the interaction between the projectile and the plasma is so weak that it 
can be considered a linear effect. In fact, one can expect that at projectile velocities 
much higher than the characteristic velocity of charged particles of the medium (e.g. the 
Fermi velocity vf), the polarizational mechanism is responsible for the main contribution 
to the energy losses. 

In order to make the calculations analytically tractable, it is usually assumed that 
the collective plasmon mode constitutes the most important channel of energy transfer 
from fast external projectiles to the plasma. If we admit the existence of a well-defined 
plasma mode with negligible damping, we can write [3] 



\\\u : (k.^ < 2 <- - 



where uif = C 2 /C , oj\ = C4/C2, C (fc), C 2 (k) and C 4 (k) are the frequency moments 
of the plasma loss function representing the zero-frequency sum rule, the f-sum rule 
and the third-frequency-moment sum rule, respectively, in completely ionized plasmas 

[H El El d . 

2. Point-like projectiles 

For a fast projectile with a point-like distribution the stopping power formula ([1]) can 
be evaluated using expression ([2]) to give [8], in a completely ionized hydrogen plasma, 

dE\ _ f Z p euj p y^ 2mv 2 (3) 



dx J cBL «»"f \ v J hujpVl + H 

h- 9 -^ 1 -^- Hp 2 SAP) d P , (4) 



67r 2 n _ 

where Z p e is the projectile charge, uj p = (Anne 2 /m) 1 ^ 2 is the plasma frequency, n is the 
electron number density, and m is the electron mass. In addition, g e i{0) is the value 
of the electron- ion radial distribution function at the origin and S e i(k) is the electron- 
ion partial static structure factor. This result implies that the well-known Bethe-Larkin 
(BL) formula [9J for the stopping power gets corrected due to the electron-ion correlation 
in the plasma (through the quantity H). In case of a one-component plasma, the 
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9ei(0) 
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1.75 


4.0 


7.4% 
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5.0 


8.7% 
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7.5 


10.5% 


4 


10.0 


11.7% 



Table 1. Modification of the Bethc-Larkin formula due to the electron- ion correlations. 
The values of g e i(0) correspond to the electron-proton plasma at 125000K, as simulated 
in Ref. [IT]. The projectiles are assumed to be protons at 1 MeV. 



existence of electron-electron correlations does not affect the fast projectile stopping 
power [ID] . 

While the correction H (or h ei (0) = g e i(0) — 1) vanishes in a weakly coupled 
plasma, when dealing with plasmas within the strong coupling regime H might become 
significant. In Table [1] the values of g e i(0) obtained by numerical simulations in Ref. 
[TT] are displayed. The simulated system was an electron-proton plasma at 125000K 
and with different values of the Brueckner parameter, r s . In the third column of Table 
[TJwe quantify, by means of the parameter 

_ {-dE/dx) cBL _ In VI + H 

(—dE/dx) bl ln^f ' 1 J 

the effect the correction H has over the stopping power asymptotic value given in ([3]), 
in this way we compare (—dE/dx) C BL to the BL formula, (—dE/dx)BL, i-e., (EI) but 
with H = 0. In addition, for the numerical estimates the projectiles were assumed to 
be protons at 1 MeV. While the results indicate that the Bethe-Larkin formula gets 
modified by some 7 -j- 11%, in a more realistic scenario one should account for other 
effects that could limit the experimental visibility of this correction. One of those factors 
is the existence of damping of the plasma mode, which would lead to the broadening 
of the plasmon peak and could modify the value of the stopping power hiding the 
correlation shift of the Langmuir mode [8]. Furthermore, the experimental techniques 
available at the present time only allow to measure the stopping power with a significant 
inaccuracy [12J. Thus, it is desirable to derive alternative methods of diagnostic where 
the correlation correction H to the stopping power could give rise to larger departures 
from its uncorrelated counterpart. 

3. Dicluster projectiles 

Let us consider then the stopping of projectiles consisting of diclusters with (in 
general, unequal) charges Zi, Z 2 whose spatial distribution (at rest) is given by 
p(r) = Zie5(r) + Z 2 e5(r — R) [2]. Under the assumption that both charge centers 
of the dicluster are randomly orientated, we can write: 

|p(k)| 2 = (Z ie ) 2 + {Z 2 ef + 2Z 1 Z 2 e 2 3o (Rk), (6) 
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Figure 1. Modification of the uncorrelated dicluster stopping power with Z\ = Z 2 due 
to the electron-ion correlation as a function of the distance between the charge centers 
of the dicluster. The same four conditions used in Table [T] are considered: r s = 1.75 
(bottom curve), r s = 2, r s = 3 and r s = 4 (top curve). 



where jo(x) is the zero-order spherical Bessel function. After substituting this charge 
distribution into Eq. ([T]) one gets [2J 

dE ( dE\ ( dE\ 

- + -3T . (7) 



dx V dx I V dx I 

\ / uncorr \ / c_. 

where the first term represents the stopping power of uncorrelated charge centers. It 
simply reduces to the same expression valid for a point-like distribution, but with the 
charge number \/ Z\ + Z\ instead of Z p . The second contribution, on the other hand, 
describes the interference pattern stemming from the correlation between the charge 
centers of the dicluster distribution [2]: 

dE\ 2e 2 nry „ f°° . ,„,^dk rkv 



2Z t Z 2 j (Rk)— cu(-lme^ (k,u))du. (8) 
Jo k Jo 



dx I TTV 2 
/ corr 

Obviously, this expression vanishes as the distance between the charge centers, R, 
increases. When R — > 0, the whole expression (JTj) reduces to the stopping power of 
a point-like distribution with the charge number [Z\ + Z 2 ). 

We can derive the fast-projectile asymptote for the correlated contribution to the 
dicluster stopping by substituting expression (j2J) into the previous formula. After 
applying the arguments presented in Ref. [8], we obtain that 

(~) * 2Z 1 Z 2 ( e -^) 2 [H(k max R)-H(k' mm R)}, (9) 
\ ax J d»"f \ v J 

\ / corr 

with TC(x) = Ci(x) — j (x), Ci(x) being the cosine integral [2], k max = 2mv/h and 
k 'min = uj p VI + H/v. 

As in the case of point-like projectiles in (jSJ), the difference between the uncorrected 
dicluster stopping power (H = 0) and the corrected one for Z\ = Z 2 is given by 



A , = In VT+H + [H(k' mm R) - H(k mm R)} 
In 5g£ + [H(k max R) - H(k min R)} 
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with k min = uj p /v. This quantity is plotted in Fig. [[J depicting its variation as a function 
of the distance between the dicluster charge centers for the thermodynamic conditions 
considered in Table [U Whereas the value of the parameter A' approaches that of A 
as R — > and R — > oo, the interference between the charge centers of the dicluster 
produces an oscillatory pattern for A' as R varies, showing that it is feasible to find 
numerical intervals for R where the dicluster stopping power exhibits a larger departure 
from the uncorrelated stopping than the one observed for the point-like stopping power 
discussed previously. For instance, in the case r s = 4 the value of A' increases from an 
11.7% to a 19%, for the optimal value of R. 

Although an extensive experimental work will be required to understand their 
practical applicability, the methods considered here might become suitable candidates 
as diagnostic tools in high energy density plasmas, where other techniques fail [12], not 
only to determine the charge number density, but even to retrieve the value of g e i{0), 
providing a method of experimental investigation of strong coupling effects in Coulomb 
systems [B]. 
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